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Abstract 

(— I ■ The infrared structure of MHV gluon amplitudes in planar limit for M = 4 super Yang- 

^ . Mills theory is considered in the next-to-leading order of PT. Explicit cancellation of the 

Oh| infrared divergencies in properly defined cross-sections is demonstrated. The remaining 

finite parts for some inclusive differential cross-sections in planar limit are calculated 
analytically. In general, contrary to the virtual corrections, they do not reveal any simple 
^ I structure. 
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O 1 Introduction 

a^ . 

• • , In recent years remarkable progress in understanding the structure of planaiEl A/" = 4 SYM 
• ^ ! (supersymmetric Yang-Mills) theory has been achieved. In the planar limit this theory seems 

^1 to be the first example of solvable non-trivial four- dimensional Quantum Field Theory. 

■ The objects which were in the spotlight starting from AdS/CFT (Anti de Sitter/ Conformal 

Field Theory) correspondence were local operators, namely the spectrum of their anomalous 
dimensions. They were calculated on the one hand side from the field theory approach [10] 
and from the other side as the energy levels of a string in classical background pi |3] revealing 
remarkable coincidence. 

Recently the quantities of interest are the so- called MHV0 scattering amplitudes. It happens 
that in the planar limit of A/" = 4 SYM theory they have a truly simple structure [H H]. It is 
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''MHV (maximally helicity violating) amplitudes are the amplitudes where all particles are treated as outgoing 
and the net helicity is equal to n — 4 where n is the number of particles. For gluon amplitudes MHV amplitudes 
are defined as the amplitudes in which all but two gluons have positive helicities 
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useful to consider the color-ordered amplitude defined through the color decomposition 



perm 



where An is the physical amplitude, are the partial color-ordered amplitudes, and aj is the 
color index of i-ih external "gluon" . 

It was found that these amplitudes reveal the iterative structure which was first established 
in two loops [5J and then confirmed at the three loop level by Bern, Dixon and Smirnov, who 
formulated the ansatz [6] for the ra-point MHV amplitudes: 



(2) 



where a = g'^Nc/4:7i and E^P vanishes as e — ^ 0, C*-'^ are some finite constants, and M^^)(/e) is 
the le- regulated one- loop n-point 0^ scalar amplitude. 

It is not surprising that the IR divergent parts of the amplitudes factorize and exponentiate 
[7]. What is less obvious is that it is also true for the finite part 



Mnie) = exp 
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where is the so-called cusp anomalous dimension [8] and Go is the second function (dependent 
on the IR regularization) which defines the IR structure of the amplitude [9] . 

According to BDS ansatz the finite part of the amplitude is defined by the cusp anomalous 
dimension and a function of kinematical parameters specified at one-loop. For a four gluon 
amplitude one has 

Fi')(0) = 2fV£ (4) 

The cusp anomalous dimension is a function of the gauge coupling, for which the four terms 
of the weak coupling expansion and two terms of the strong coupling expansion [21 [3] are 
known. Integrability from the both sides of the AdS/CFT correspondence leads to all-order 
integral equation [llj solution of which being expanded in the coupling reproduces both the 
series [12]. 

While for n = 4, 5 the BDS ansatz goes through all checks, namely the amplitudes were 
calculated up to four loops for four gluons and up to two loops for five gluons [13J. However, 
starting from n = 6 it fails. The first indication of the problem was strong coupling calculation 
in the limit n —>■ oo where discrepancy with the BDS formula was found. The second 
indication came from the amplitude/Wilson loop duality [151 lini El HE] namely from the 
comparison of the hexagonal light-like Wilson loop and finite part to the BDS ansatz for the 
six-gluon amplitude. It was found that the two expressions have difference in a non-trivial 
function of the three (dual) conformally invariant variables [I9] . The third indication appeared 
in the paper [20] where the analytical structure of the BDS ansatz was analyzed and starting 
from n = 6 the Regge limit factorization of the amplitude in some physical regions fails. Finally 
it was shown by explicit two-loop calculation [21] that the BDS ansatz is not true and it needs 
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to be modified by some unknown finite function, which is an open and intriguing problem. 
From two-loop calculation for the six-point amplitude [21J and hexagonal light-like Wilson loop 
[22] it was shown that the gluon amplitude/ Wilson loop duality is still valid. 

While all the UV divergences in A/" = 4 SYM are absent in scattering amplitudes the IR 
ones remain and are supposed to be canceled in a properly defined quantities. Regularized 
expressions act like some kind of scaffolding which has to be removed to obtain eventual phys- 
ical observable. It is these quantities that are the aim of our calculation. And though the 
Kinoshita-Lee-Nauenberg [23] theorem in principle tell us how to construct such quantities, 
explicit realization of this procedure is not simple and one can think of various possibilities. In 
particular, one could consider the so-called energy flow functions defined in terms of the energy- 
momentum tensor correlators introduced in [21] and considered at weak coupling regime in [25] 
and recently at strong coupling regime in [26]. From our side we concentrated on inclusive 
cross-sections in the hope that they reveal some factorization properties discovered in the reg- 
ularized amplitudes. Similar questions were discussed in [23, where the inclusive cross-section 
like the IR safe observables in A/'=4 SYM were constructed. 

2 The Infrared Safe Observables 

To perform the procedure of cancellation of the IR divergences one should have in mind that 
in conformal theory all the masses are zero and one has additional collinear divergences which 
need special care. In these work we employ the method developed in QCD parton model 
[281 [291 l30l [3T[ 132] . It includes two main ingredients in the cancellation of infrared divergencies 
coming from the loops: emission of additional soft real quanta and redefinition of the asymptotic 
states resulting in the splitting terms governed by the kernels of DGLAP equations. The latter 
ones take care of the collinear divergences. 

Typical observables in QCD parton model calculations are inclusive jet cross-sections, where 
the total energy of scattered partons is not fixed since they are considered to be parts of the 
scattered hadrons. In [21] the algorithm of extracting divergences was developed which allows 
one to cancel divergences and apply numerical methods for calculation of a finite part. In our 
paper we choose as our observables the inclusive cross-sections with fixed initial energy and get 
analytical expression for the finite part of the differential cross-section. We do not assume any 
confinement and consider the scattering of the single parton based coherent stated being the 
asymptotic states of conformal field theory. 

When the number of particles increases one has to specify the measurable quantity and to 
distinguish the particle(s) in the final state. If one wants to construct the finite quantity it 
is not sufficient to consider the process with the fixed number of final particles. One has to 
include processes with emission of additional soft and collinear massless states, i.e to consider 
the inclusive cross-section. 

One possibility is to introduce the energy and angular resolution of the detector and to 
cut the phase space so that the soft quanta with total energy below the threshold as well as 
all the particles within the given solid angle are included. This procedure works well in QED 

^squared perturbative amplitudes used in our calculation have been summed over colors, so in this sense they 
are colorless and there are no contradiction with statements that cancellation of IR divergences occurs only for 
colorless objects. 
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but introduces explicit dependence on the energy and angular cut off, thus violating conformal 
invariance. 

We adopt here the other attitude and do not introduce any cut off in the phase space but 
rather consider the inclusive cross-section with emission of all possible particles allowed by 
kinematics. Then one has to specify which particle is detected. For instance, one can measure 
the scattering of a given particle on a given angle integrating over all the other particles. As 
it will be clear later in this case one still cannot avoid introducing some scale related to the 
definition of the asymptotic states of a theory. 



3 Calculation of Inclusive Cross-sections in A/^ = 4 SYM 
theory 

Our aim is to evaluate the NLO correction to the inclusive differential polarized cross section in 
the weak coupling limit in planar A/" = 4 SYM in analytical form and to trace the cancellation 
of the IR divergences. 

We start with the 2^2 MHV scattering amplitude with two incoming positively polarized 
gluons and two outgoing positively polarized gluons and consider the differential cross-section 
da2~,2{g^g^ g^g^)/d^ as a function of the scattering solid angle. The total cross-section 
is divergent at zero angle. Treating all the particles as outgoing this amphtude is denoted as 
( , ++) MHV amplitude. At tree level the cross-section is given by 



\Mr>VS2. (5) 



where J is flux factor, in our case J = s, s is standard Mandelstam variable, (i02 is the 
phase volume of two-particle process (we use dimensional regularizatiori§, D = A — 2e) and 
Sn {n = 2 in this case) is the so-called measurement function which specifies what is really 
detected, cil^is = (i0i3(icos(6'i3)lll, 6^13 is the scattering angle of the particles with momenta 
P3 with respect to in the center of mass frame. The matrix element is obtained from the 
color-ordered amplitudes via summation [H H] 

|_^{tree)(-++)|2 ^ ^4^2(^2 fl2 ^ (g) 

where in all expressions we take Sij = {pi +pj)'^. 

Within dimensional regularization(reduction) the cross-section looks like 
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(7) 



where s, t, u are the Mandelstam variables, E is the total energy in the center of mass frame, 
c = cos6'i3, yU and e are the parameters of the dimensional regularization. The next step is to 
calculate the NLO corrections. 



^ namely - FDH version of dimentional reduction, see for details [33] 

^if to be more accurate in dimensional regularization(reduction) we have (Kl\r^ 
bi3sin{(j)i3)~^''dcos{9i3)sin{0i3)~^'' . 
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3.1 Virtual part 



To get the one-loop contribution to the differential cross section we use already known one loop 
contribution for the color-ordered amphtude 



(l~loop) I 1 (tree) 



1 



where I'f^i^s^t) is the scalar box diagram 



r(l-/oop) 



2 rri+e)rri-e) 



--sth 
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r l-2e 



-T -r+ — r +:.iogM — 



-t 



Then it is straightforward to obtain the one-loop contribution to the cross-section in the 
planar limit 



(— ++) 



-^■J / v'lrt 

+— ^logi 



1 + c, 



2e 



a 



16 ?> + & 
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16(3 + c^)7r^ 
3(1 -c2)2 
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,2^2 



log I 



1 -c. 



(1-C 



,2^2 
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It should be stressed that due to conformal invariance of A/" = 4 SYM theory at quantum level 
there are no UV divergences in (JHl) and all divergences have the IR or coUinear nature. They 
have to cancel in the properly defined observables. 



3.2 Real emission 

The next step is the calculation of the amplitude with three outgoing particles. Here we have 
to define which is the process that we are interested in. There are several possibilities. 

1. Three gluons with positive helicities: g^g^ —>■ g^g^g^ ■ This is the MHV amplitude; 

2. Two gluons with positive helicities and the third one with negative helicity: g^g^ — > 

1^ This is the anti-MHV amplitude; 

3. One of three final particles is the gluon with positive helicity and the rest is the quark- 
antiquark paiil§: g'^g^ g^Q Q^ oi g^g^ g^Q^Q - This is an anti-MHV amplitude; 

4. One of three final particles is the gluon with positive helicity and the rest are two scalars: 
g~^g~^ g~^AA. This is an anti-MHV amplitude. 

If one fixes one gluon with positive helicity scattered at angle 6 and sum over all the other 
particles then all the processes mentioned above contribute. In case when one fixes two gluons 
with positive helicity and looks for the rest, only the first two options are allowed. 

^There is also a g^g^ g^g~ 9^ helicity configuration, but the partial amplitudes for them are equal. We 
will use ( h H — ) notation for both of them. 

^The A/" = 4 supermultiplet consists of a gluon g, 4 fermions ("quarks") and 6 real scalars A"^-^. A and 
B are SU (4)ij indices, A is an antisymmetric tensor. It is implied that all squared amplitudes with quarks and 
scalars are summed over these indices. 



5 



The cross-section of these processes can be written as 



da' 



2^3 
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1 
J 



(tree) |2 



5, 



(9) 



where dcj)^ is the three-particle phase volume and is the measurement function, which con- 
straints the phase space and defines the particular observable. 

We omit the details of the calculation and for the lack of space present here only the divergent 
parts of the calculated objects. We leave all the full answers for a separate publication and 
present the finite part only in the simplest case below. 

When there are three identical particles in the final state one has to define which ones are 
detected. In case of one detectable particle one can choose the fastest one, in case of two - the 
two fastest. Selection of detectable particles can be achieved restricting the momentum. Thus 
choosing > E/3 we fix the gluon with momentum as the fastest particle. When the final 
particles are not identical this problem does not appear and the phase space is not restricted. 
In what follows we restrict the moment of the gluon by a universal value p^ > {1 — 6)/2E and 
keep the value of 6 arbitrary. The case of identical particles then corresponds to 5 = 1/3 and 
the case of nonidentical particles to 5 = 1. One can show that in this case the requirements of 
stability of observable with respect to emission of soft and collinear quanta [22] are satisfied. 
We show below that IR and collinear divergences cancel in observables for an arbitrary value 
of (5. 



da 



1. Real Emission (MHV) 

(—+++) „,2 at2 



dn 




Notice the singularity when 6^1. 
2. Real Emission (anti-MHV) 

-++-) 




Real 




12(c2 + 3)log5 64(12c2 + 17) 



+ 



3(1 



;(5 + 3c^)6^ - (c^ + 19)5 + 2(5c2 + 43) + 



[l-d^y V3 

8(c2-6c + 21), ,l + 5-(l-5)c. 



'2(3c2-24c + 85) , A 



;i-c)(l + c)3 



log( 



+ 



(1-C)(1 + C)3 

32(2- 



log( 

c) 



32(c2-4c + 7) 



(l + c)3(l 
64(1 - c) 



c)(l 



:i + c)3(l + 5 
3. Fermions 



c(l-5))2 3(l + c)3(l + 5-c(l-5))^ 



c{l-6)) 



+ (c ^ -c) 



'11^ 



+ Finite part 




(12) 



+ Finite part 
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4. Scalars 



' da- 



2->3 



(— +AA) 



'13 / Real 

3(5 - c) 



^2 



:i + c) 



log I 



1-c, 3(5 + c), ,1 + c' 
- 7^ ITT log(^^) 



2(10 + 7c2 
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3.3 Splitting 

Taking into account emission of additional soft quanta allows one to cancel the IR divergences 
(double poles in e) but leaves the single poles originating from collinear ones. Indeed, in case 
of massless particles the asymptotic states (both the initial and final ones) are not well defined 
since a massless quanta can split into two parallel ones indistinguishable from the original. To 
take this into account one introduces the notion of distribution of the initial particle (gluon) 
with respect to the fraction of the carried momentum z: g{z). Then the initial distribution 
corresponds to g{z) = 6{1 — z), and the emission of a gluon leads to a splitting: the gluon 
carries the fraction of momentum equal z, while the collinear gluon - {1 — z). The probability 
of this event is given by the so-called splitting functions Pgg{z). In case of a gluon in a final 
state this corresponds to the fragmentation of the gluon into pair of gluons or pair of quarks 
or scalars. 

Additional contributions from collinear particles in initial or final states to inclusive cross- 
sections have the form, respectively 

,2 



dcr. 



spl,init 
2-»2 




dzPgg{z) J2 

«J=1,2; i^j 



da2^2{zPi,Pj,P3,P4)S2 



spl,init I 



*2iPl,P2,P3,Pi) i dz J2 Pgl{z)Sf'^'''{ 



(14) 



(15) 



l=g,q,A 



where the scale Q'j, sometimes called the factorization scale, belongs to the definition of the 
coherent asymptotic state and restricts the value of transverse momenta. The dependence of 
parton distribution on Q'j is governed by the DGLAP equation. The splitting function Pij 
for each helicity configuration can be obtained as a collinear limit of the corresponding partial 
amplitude (see for example |3B] fore more details). 

Taking into account the splitting of initial states and the fragmentation of the final states 
we get the following contribution to the inclusive cross sections 

1. The initial and final splitting for the MHV amplitude. 

(— + 



da- 



2^3 



dn 



13 / InSplit 

8(c2 + 3) 

fl-c2)2 




' da- 



2^3 



dn 



13 



(—+++) 



FnSplit 




1 4(c2 + 3) , 1 
log- 



e(l-c2)2 « 5 



(17) 
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2. The initial and final splitting for the anti-MHV amplitude 
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(1 - C2)2 
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+ 



(1 - C2)2 

8(c2-6c + 21) l + 5-c{l-5) 
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^log- 
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(l + c)3(l-c)(l + (5-c(l-(5)) 
64(1 - c) 
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+ (c -c) 



+ Finite part | , 



' da^ 



2^3 



(—++-) 



a'N^fi'Y fx'Va ri4(c^ + 3) 



dn 



13 / FnSplit 



S )\Qj) TT \e (l-c2)^ 



3. The initial splitting for the quark final states (5—1) 



d(J2^3 

, d^i3 



( +qq) 



+ 



InSplit 

2(3 -c)2 



E'^ 
1 - 



« J 16 
TT 1 e 



log5-5(^52-^(5+3) 



(79 + 25c2 



(19) 



, c, 2(3 + c)2 



(1-C)(1 + C)3 2 ^ ' (l-c)3(l + c) 

4. The initial splitting for the scalar final states (5 = 1) 



log( 



3(1-C2)2 
1 + C 



(20) 
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' da^ 



2^3 
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2(10 



3(5 -c) 1-c, 
■log(- 



(1 - C2)2 

+ Finite part > . 



(21) 



3(5 + c) 1 + c 
(1 + cf ^^''^ 2 ' (1 - c)3 2 ^ 

4 IR safe observables in A/" = 4 SYM 

In the NLO there are two sets of amplitudes, namely the MHV and anti-MHV amplitudes, 
which contribute to the observables. The leading order 4-gluon amplitude is both MHV and 
anti-MHV and we split it into two parts. Then one can construct three types of infrared-safe 
quantities in the NLO of perturbation theory, namely 

• pure gluonic MHV amplitude 



A 



MHV 



1 I da2^2 

2 ~ 



dQ 



(— ++) 
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13 / virt 
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pure gluonic anti-MHV amplitude 



B 



anti-MHV 



1 / d(72_^2 
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dn 



(--++) 
I + 
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13 / InSplit 
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; (23) 
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anti-MHV amplitude with fermions or scalars fropm the full Af = 4 supermultiplet 



C 



Matter 



'da- 



( + , qq+Kk) 
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da 



2^3 



13 / Real 



dn 



13 
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InSplit 



(24) 



We would like to stress once more that in each expression fl22f23f24p all IR divergences cancel 
for arbitrary S and only the finite part is left. 

Defining now the physical condition for the observation we get several infrared-safe inclusive 
cross-sections 



Registration of two fastest gluons of positive hehcity 



A 



MHV 



5=1/3 



+ B 



anti-MHV 



5=1' 



Registration of one fastest gluon of positive helicity 



A 



MHV 



(5=1/3 



+ B 



anti-MHV 



5=1/3 



+ c 



Matter 



5=1' 



(25) 



(26) 



Anti-MHV cross-section 



B 



anti-MHV 



5=1 



c 



Matter 



5=1 



(27) 



Relative simplicity of the virtual contribution ([H]) which does not contain any special func- 
tions but logs suggests similar structure of the real part. However this is not the case. While the 
singular terms are simple enough and cancel completely the finite parts are usually cumbersome 
and contain polylogarithms. The only expression where they cancel corresponds to 5 = 1 case 
which is possible only for the last set of observables, namely for the anti-MHV cross-section 
Choosing the factorization scale to be Q/ = we get 




AntiMHV 



3 + c2 
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,2^2 
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3(l + c)3(l -c)^ 
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One can see that even this expression does not repeat the Born amplitude and does not have any 
simple structure. Note, that the finite answer depends on the factorization scale. This depen- 
dence comes from the asymptotic states which violate conformal invariance of the Lagrangian. 
This dependence seems to be unavoidable and reflects the act of measurement. Construction 
of observables which do not contain any external scale remains an open question. 
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5 Discussion 



To solve the model might have different meaning. Calculation of divergences and understanding 
of their structure is very useful but surely not enough. The knowledge of the S-matrix would 
be the final goal though the definition of the S-matrix in conformal theory is a problem. Even 
in the absence of the UV divergences there are severe IR problems and matrix elements do not 
exist after removal of regularization. The experience of QCD, which is very similar to A/" = 4 
SYM theory from the point of view of the IR problems, tells us that in inclusive cross-sections 
the IR divergences cancel and one has finite physical observables. However, one either has to 
redefine the asymptotic states or consider the scattering of "hadrons". In both the cases one 
has to introduce some parton distributions which are the functions of a fraction of momenta 
and, in higher orders, of momenta transfer. This leads to appearance of a factorization scale 
which breaks conformal invariance. This means that we do not have meaningful observables in 
a pure conformal theory. The finite observables of the type considered here are the inclusive 
cross-sections 



which besides the kinematical variables contain the dependence on the factorization scale. 

Remarkable factorization properties of the MHV amplitudes accumulated in the BDS ansatz 
(with the so far unknown modification) suggest the way of receiving "exact" results for the 
amplitudes on shell. However, as we have already mentioned, it is the finite part that we 
are really for. Unfortunately, our calculation has demonstrated that the simple structure of 
the amplitudes governed by the cusp anomalous dimension has been totally washed out by 
complexity of the real emission matrix elements integrated over the phase space. This means 
that either J\f = 4 SYM theory does not allow such a simple factorizable solution or that we 
considered the unappropriate observables. Alternatively one may relay on simplification in 
particular kinematic regime like the Regge limit where the exponentiation is expected. 

There is an interesting duality between the MHV amplitudes and the Wilson loop, between 
the weak and the strong coupling regime [ISl IISI El [IE]- Probably it would be possible using 
the AdS / CFT correspondence to construct the IR safe observables at the strong coupling limit 
(similarly to what we did here) and to shed some light on the "true" calculable objects in 
conformal theories. 
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